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Abstract 

We analyze the partition function of the dimer model on an A4 x J\f triangular lattice wrapped 
on torus obtained by Fendley, Moessner and Sondhi [Phys. Rev. B 66, 214513 (2002)]. From 
a finite-size analysis we have found that the dimer model on such a lattice can be described by 
conformal field theory having central charge c = — 2. The shift exponent for the specific heat is 
found to depend on the parity of the number of lattice sites J\f along a given lattice axis: e.g., for 
odd J\f we obtain the shift exponent A = 1, while for even M it is infinite (A = oo). In the former 
case, therefore, the finite-size specific-heat pseudocritical point is size dependent, while in the latter 
case, it coincides with the critical point of the thermodynamic limit. 
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I. INTRODUCTION 



In experiments and in numerical studies of critical phenomena, it is essential to take 
into account finite-size effects in order to extract correct infinite-volume predictions from 
the data. Therefore, in recent decades there have been many investigations of finite-size 
scaling, finite-size corrections, and boundary effects for critical model systems. In the quest 
to improve our understanding of realistic systems of finite extent, two-dimensional models 
play a crucial role in statistical mechanics as they have long served as a testing ground to 
explore the general ideas of finite-size scalingunder controlled conditions. Very few of them 
have been solved exactly, the Ising model (l-9| and the dimer model 10- 15| being the most 
prominent examples. 

The dimer model is a two-particle system. The main difference between it and one-particle 
systems such as Ising, Heisenberg, or Potts models etc., is that occupation of a given lattice 
site ensures that at least one of its nearest-neighbor sites is also occupied. It is well known 
that, due to this non- locality, the critical behavior exhibited by the dimer model can depend 
on the lattice structure and shape (square, triangle, honeycomb, etc). Previous studies have 
shown that finite-size corrections in the free energy can exhibit a strong dependence upon 
the parity of the lattice, and this has provoked controversial conclusions about the value 
of the central charge from c = — 2 to c = 1. One can expect that such unusual finite-size 
behavior should also hold for the dimer model on other lattices and here we investigate the 
triangular lattice in particular. 

We are particularly interested in the finite-size scaling behavior of the specific-heat pseu- 
docritical point. In finite systems the counterparts of the singularities which mark higher- 
order phase transitions in the thermodynamic limit are smooth peaks the shapes of which 
depend on the critical exponents. In particular, let C (t, L) be the specific heat at a re- 
duced temperature given by t for a system of linear extent characterized by L. In the 
infinite- volume limit, C(t, oo) diverges at the critical point t — t c — 0. In finite volume, 
the analog to the divergence is a finite peak the shape of which is characterized by (i) its 
position tp S eudo (h) its height C(t pseu d , L) and (iii) its value at the infinite- volume critical 
point C(0, L). In particular, the position of the specific heat peak, t pscu do, is a pseudocritical 
point which approaches t c = as L~ x , where A is called the shift exponent. In most models 
exhibiting higher-order phase transitions, the shift exponent coincides with the inverse of 
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the correlation-length critical exponent 1/u, but this is not a direct conclusion of FSS theory 
is not always true. 

For example, for the Ising model in two dimensions, Ferdinand and Fisher determined 
that behavior of the specific-heat pseudocritical point matches that of the correlation length 
with A = l/z/=l[l[. However, Ising models defined on two-dimensional lattices with 
other topologies have shift exponents which differ from the inverse correlation length critical 
exponent (see Ref.Q and references therein). This is despite the fact that the critical 
properties on such lattices are the same as for the torus in the thermodynamic limit. A 
question we wish to address here is the corresponding status of the shift exponent in the 
dimer model. 

In contrast to spin models, the critical behavior of dimer models are strongly influenced 
by the structure of the underlying lattice. For example the square lattice dimer model is 
critical with algebraic decay of correlators Jl6l. while the dimer model on the anisotropic 



18| . exhibits 



honeycomb lattice, which is equivalent to five- vertex model on the square lattice 
a KDP-type singularity and the dimer model on the Fisher-type lattice exhibits Ising-type 
transitions jl9|. Thus, it appears that the dimer model itself has not a single critical behavior, 
but several critical behaviors associated with different classes of universality. 

It has been shown explicitly 20| that the free energy per site for the dimer model on 



the square lattice is insensitive to the precise form of the boundary conditions in the limit 
of a large lattice. This is in contrast to its finite-size counterpart, for which sensitivity to 
boundary conditions is notable feature, in particular to the parity of the number of lattice 



sites along a given lattice axis 
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2l| . Similar statements hold for the dimer model on the 



honeycomb and triangular lattices. 
Very recently, it has been shown 



12j that the finite-size corrections of the dimer model 



on planar oo x M square lattices also depend crucially on the parity of M and the boundary 
conditions and such unusual finite-size behavior can be fully explained in the framework of 
the c = — 2 logarithmic conformal field theory. 

Our objective in this paper is to study the finite size properties of dimer model on the 



plane triangular lattice using the same techniques developed in papers [4|, The paper 
is organized as follows. In Sec. II we introduce the dimer model on the triangular lattice 
with periodic boundary conditions. In Sec. Ill we discuss the finite-size corrections for an 
infinitely long cylinder of circumference M and find that the dimer model on the triangular 



3 



> 



Figure 1: The Ai x M triangular lattice with dimer weights Zh in the horizontal direction, z v in the 
vertical direction, and t in the diagonal direction. 

lattice can be described by conformal field theory with central charge c = —2. In Sec. IV 
we investigate the properties of the specific heat near the critical point and find that the 
specific-heat shift exponent A depends on the parity of the number of lattice sites along the 
lattice axis N '. For odd M we obtain for the shift exponent A = 1, while for even M we find 
that the shift exponent A is infinity (A = oo). Our results are summarized and discussed in 
Sec. V. 

II. PARTITION FUNCTION 

In the present work, we consider the dimer model on M. x M triangular lattice under 
periodic boundary conditions. The partition function is given by 



where summation is taken over all dimer covering configurations, where Zh, z v and t are, 
respectively, dimer weights in the horizontal, vertical and diagonal directions, and where n^, 
n v and n t are, respectively, the number of horizontal, vertical and diagonal dimers (FigJTJ). 
The dimer model on the triangular lattice undergoes phase transition at the point t — t c — 
(and likewise for Zh and z v ), where the partition function is singular. Thus the general 
triangular lattice model is critical in the square lattice limit. The dimer weight t plays a 
role similar to the reduced temperature in the Ising model. In what follows, we will set 
Zh — z v — 1- 

An explicit expression for the partition function of the dimer model on an ]V[ x M 




(1) 
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triangular lattice wrap ped on torus has been obtained by Fendley, Moessner and Sondhi 22 j 
and can be written as 



ZmM*) = 2 [Gofi(t,M,AT) + G 0A/2 (t,M,AT) + G 1/2fi (t, M,Af) + G 1/2jl/2 (t,M,N)] 



(2) 



where 



M/2-1M-1 



Gl^MM)= 1 4 



sin 



m=0 n=0 

2 2vr(m + f$) 



2n{n + a) 
77 



+ sm 



i'2ir(n + a) 2ir(m + 8) 

t cos 



M ' V M ' M 

for even Ai. The notation a = corresponds to periodic boundary conditions for the 
underlying free fermion in the TV-direction while a = \ represents anti-periodic boundary 
conditions. The boundary conditions in the A^-direction are similarly controlled by the 
parameter /3. 

Since the total number of sites must be even if the lattice is to be completely covered by 
dimers, we will consider two cases, namely even-even (ee) case when Ai = 2M and M = 2N, 
and even-odd (eo) case when Ai = 2M and M = 2N + 1. Note that due to the symmetry of 
the lattice the odd-even case (oe) (Ai = 2M + 1, M = 2N) can be obtained from even-odd 
case by simple transformation £ — > l/£, where £ = AijM is a aspect ratio. 

a). Dimers on 2M x 27V lattices 

In the even-even case where (Ai,Af) = (2M,2N), the second product in ([3]) may be 

2N-X 

compactly written as [ F(n + a, m + 0) where the function F(x, y) is given by 



n=0 

F(x,y) =4 



• 2 ™ , • 2 W , 2 2 /vrx ny 
Sm lV +Sin M +t C ° S [JT + M 



(4) 



Splitting this product into two parts, 

2JV-1 N-l 2N-1 

H F(n + a,y)=l[F(n + a,y)x + (5) 

n=0 n=0 l=N 

shifting the index in the second part from I to n = I — N, and noting the translation 
symmetry F(N + x) = F(x), this may be expressed as 



2AT-1 



'N-l 



[ F(n + a,y) = [ [F(n + a,y) 



(6) 



n=0 



,n=0 



Defining the partition function with twisted boundary conditions 

M-l N-l 



Z 2 a , (t,M,N) = J] II 4 



m=0 n=0 



sin' 



ir(n + a) 
TV 



. 2 7r(m + /3) 2 2 /7r(ra + a) 7r(m + /3) 
+ sin h t cos 



Af 



iV 



M 



(7) 



one has 



Ga, p (t, 2M, 2N) = Zlp(t, M, N). (8) 
The even-even partition function given by Eq.(j2J) and Eq.([3]) can now be written in the form 
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Z 2M Mt) = \ [2%, (t t M, N) + Z* tl/2 (t, M, N) + Z? /2i0 (t, M, N) + Z* /2>1/2 (t, M, N)] . (9) 

Note that Eq.([71) at t — coincides with the corresponding expressions for the square lattice 
for which a general theory about its asymptotic expansion has been given in Ref . |^ | . 

Note also that for all (a, (3) ^ (0,0), the partition function Z a ^(t, M, N) is even with 
respect to its argument t. Hence, near the critical point (t = 0) we have 

Z a , p (t,M,N) = Z aiP (0,M,N) + t ^Z" at p(0,M,N) + ... for (a, (5) ^ (0, 0). (10) 

The only exception is the point where both a and /3 are equal to zero. This case has to be 
treated separately since at this point (t = 0) the partition function Z ^(t, M, N) vanishes. 
As a result, we have 

t 3 



Z 0fi (t, M, N) = tZ Ot0 (0, M, N) + -Z'" 0i0 (0, M, N) 
At the critical point t = we have 



for (a, /3) = (0,0). (11) 



Z 0>0 (Q,M,N) = 



JV-l 



Z Qi/3 (0,M,A) = ' 



n=0 



2 sinh 



, 7r(n + a) . 



(12) 

for («,/3)^(0,0), (13) 



where u(k) = arcsinh (sin k). In the derivation of the Eq. (fT3|) we have used the identity 23 1 



M-l 



m=0 



sinh 2 o; + sin 2 (^±£)) 



4 | sinh (Mw + i%/3)\ 2 . 



(14) 
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Taking the derivative of Eq.(j71) with respect to variable t and then considering limit t — > 
we obtain 

AT-1 

Z' 0fi {0,M,N) = 2M J]2|sinh(Mu;(^))|, (15) 

n=l 

Z' aiP (0,M,N) = for (a,/?) ^(0,0). (16) 
b). Dimers on 2M x (2iV + 1) lattices 



In Ref. 



15] , it has been shown that in the even-odd case, the partition function given by 



Eq-d2D and Eq.([2]) can be written as 

Z 2M ,2N+i(t) = Z 0)0 (i, M, 2N + 1) + Z 0tl/2 {t, M, 2N + 1). (17) 

Note that Eqs. and (j!7p in the case t — coincide with the corresponding expressions 
for the square lattice (see Ref. Jll|). 

Thus we can see that partition function for the dimer model on triangular lattice under 
periodic boundary conditions can be expressed in terms of the only one subject, namely, 
Z a ^(t, M, N) with (a, 0) = (0, 0), (0, 1/2), (1/2, 0) and (1/2, 1/2). 

III. DIMER ON THE INFINITELY LONG STRIP 

Conformal invariance of the model in the continuum scaling limit would dictate that 
the asymptotic finite-size scaling behavior of the critical free energy of an infinitely long 
two-dimensional strip of finite width has the form 

/ = + ^ + ^5 + ..., (18) 

where /bulk is the bulk free energy, / sur f is a surface free energy and A is constant. Unlike the 
free energy densities and / sur f, the constant A is universal. The value of A is related 
to the central charge c and the highest conformal weight A of the underlying conformal 
theory, and depends on the boundary conditions in the transversal direction. These two 



2J-|26|, 



dependencies combine into a function of the effective central charge c e ff = c — 24A 

7T / C \ 

A = -7^c e// = 7T^A-—J on a strip, (19) 

7T / C \ 

A = ——c e ff = Air ( A — — J on a cylinder. (20) 



Let us now consider the dimer model on the infinitely long strip of width M under periodic 
boundary conditions. 

Considering the logarithm of the partition function given by Eq. tTTBl . we note that it can 
be transformed as 

JV-l N-X 



In Z Q>/3 (0, M, N) = M <» (^P) + E ln 



(21) 



n=0 n=0 

The second sum here vanishes in the formal limit M — > oo. The asymptotic expansion of 
the first sum can be found with the help of the Euler-Maclaurin summation formula 

N — 1 r, oo / 2 \ P \ T3a 

M^(^)=ll^-^-^Y,(^) <^£*|, (22) 

where J q lu(x) dx = 27, 7 = 0.915965... is Catalan's constant and B° are so-called Bernoulli 
polynomials and S = MN. We have also used the symmetry property, u(k) = uj(tt — k), of 
the lattice dispersion relation u(k) and its Taylor expansion 

00 , 

"<*> = (23) 

where A = 1, A 2 = —2/3, A 4 = 4, etc. 

Thus one can easily write down all the terms of the exact asymptotic expansion for the 
F a A N ) = ~ limM^oo jj In Z a>/3 {M, N) 

F a JN) = - lim -\nZ a MN) = -^N + 2y(-) 2P+1 ^-^. (24) 

From F a ^(N), we can obtain the asymptotic expansion of the free energy per bond of an 
infinitely long cylinder of circumference M . Since the expression for the partition function 
is different for even M and odd M, we will consider these two cases separately. For even M 
(N = 2N), we have 

f = " IMn 111 Z ™*»M = " J3L JMn ln z ^" iM - N) = a^wW. (25) 

and for odd Af 
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mToo 2M(2N + 1) " mToc 2M(2~N + 1) 



f = ~ .i im ^777777 I ^ ln^2M,2JV+i(0) = - .ton ollif , OAr , , , In Z , 1/2 (M, 2N + V 



F ,i /2 (2JV+1). (26) 



2(2iV+ 1) 



From Eg. (1251) using Eq. (l24]) one can easily obtain that for even M the asymptotic expan- 
sion of the free energy is given by 

\ 2p+2 ^ R l/2 



p=0 x ' 



^2p J - , 2p+ 2 

(2p)\ 2p + 2 



= /buik-^ + ... (forAT=2iV), (27) 
while for odd Af from Eqs. (|26|) and (124)) one can obtain 

/ - /bulk + ^2^\m) (2p)!2p + 2 

p=0 

= /bulk + "... (for A/" = 2iV + 1). (28) 

The bulk free energy /bulk = — ^ is the same for M even and odd cases. Thus we find that 
the finite-size corrections in a crucial way depend on the parity of M. In particular it means 
that due to the certain non-local features present in the dimer model, a change of parity of 
M induces a change of boundary condition. The similar situation also happen in the dimer 



model on the square lattice, see 12[, where a detailed analysis of boundary conditions and 
parity dependence effects has been carried out in this context. 

Since the effective central charge merely determines some combination of c and A, one 
cannot obtain the values of both without some assumption about one of them. This as- 
sumption can be a posteriori justified if the conformal description obtained from it is fully 
consistent. Surprisingly, there are two consistent values of c that can be used to describe 
the dimer model, namely c = —2 and c — 1. For example for the dimer model on an 
infinitely long cylinder of even circumference M one can obtained from Eqs.( fl8l) . ( 1201) and 
( 1271) that the central charge c and the highest conformal weight A can take values c = 1 
and A = or c = —2 and A = —1/8. For the dimer model on an infinitely long cylinder of 
odd circumference M one can obtained from Eqs. ( I18[) . (120]) and ( 128]) that the central charge 
c and the highest conformal weight A can take values c = 1 and A = 1/16. It turns out 
in this case that another consistent conformal description exists, with c = —2 and A = 



|12| . In particular, it has been shown (for more details see [12]) that although the dimer 
model is originally defined on a cylinder with odd circumference A/", it shows the finite-size 
corrections expected on a strip and must really be viewed as a model on a strip. 

Thus from the finite size analyzes we can see that two conformal field theories with the 
central charges c = 1 and c = —2 can be used to described the dimer model on the triangular 



lattice. But since the general triangular lattice model is critical in the square lattice limit 



12|, 



we come 



and the dimer model on the square lattice belongs to c = —2 universality class 
to the conclusion that the dimer model on the triangular lattice can also be described by 
conformal field theory having central charge c = —2. 



IV. SPECIFIC HEAT NEAR THE CRITICAL POINT 

Let us now consider the behavior of the specific heat near the critical point. The specific 
heat C(t, Ai,J\f) of the dimer model on KA x M triangular lattice is defined as 

C(t,M,M) = --^f(t,M,Af), (29) 

where f(t,A4,Af) is free energy of the system 

f(t,M,Af) = -~hxZ M Mt)> (3°) 

and where S = M.M is the lattice area. 

The pseudocritical point t pscu do is the value of the temperature at which the specific 
heat has its maximum for finite M. x M lattice. One can determine this quantity as the 
point where the derivative of C(t, Ai,J\f) vanishes. The pseudocritical point approaches the 
critical point t c = as L — > oo in a manner dictated by the shift exponent A, 

l^pseudo — t c \ ~ L A . (31) 

where L = \/S is the characteristic size of the system. The coincidence of A with where 
v is the correlation lengths exponent, is common to most models, but it is not a direct 
consequence of finite-size scaling and is not always true. 

Since the expression for the partition function is different for even M and odd M, we 
consider these two cases separately. 

Let us start with case of odd M (N = 2N + 1). Expanding the expression (129]) about 
the critical point t — with the help of Eqs.f fTOj) . f TTTT) and (fl7|) yields 



t 2 

C{t,M,M) = C(0,M,Af) + tC {l) (0,M,Af) + — C (2) (0, M,N) + 0(t 3 ), (32) 
where C{0,M,Af) is the critical specific heat, and C (n) (0, M, Af) = —C(t,M,Af)\ t =o- We 
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have 

r(2) 



SC V> M >*n - Z , 1/2 (0,M/2,A0 {z , 1/2 (0,M/2^)) ' (33) 

qr(1)((] M .rs = Z<$(0,M/2,Af) o Z^(0,M/2,Af)Z^ /2 (0,M/2,Af) 
l ' ' J Z 0)1/2 (0,^/2,^) Z 2 1/2 (0,M/2,AT) 



Z 0tl/2 (0,M/2,Af) 



(34) 



9 r<%n M An g^(g>VVVO , „ / zg)(0,M/2,JSQ \ Z$ /2 (0,M/2,M) 
1 ' ' J "' Z 0il/2 (0,M/2,M) + \Z 0tl/2 (0,M/2,M)J Z , 1/2 (0,M/2,Af) 

Z$ /2 (0,M/2,M)\ 2 Z$(0,M/2,M)Z$(0,M/2,M) 



3 



Z 0il/2 (0,M/2,Af) / Z* 1/2 (0,M/2,M) 



I g^M^L ^i (35) 

From Eg. (|32|) . the first derivative of the specific heat on a finite lattice near the infinite 
volume critical point can be found, and it seen to vanish when 

To find the exponent A one need the finite-size corrections to C^'(0,/A,M) and 
C^ 2) (0,M,Af). The exact asymptotic expansions of C w (0,M,Af) and C (2) (0, .M,TV) can 



be found along the same lines as in Ref . |27| and leading finite size behavior is 

C m (0,M,Af) ~ v 7 ^, (37) 
C {2) (0,M,N) ~ 5. (38) 

Expressions Eq. (l37p and (I3"8"j) now gives the FSS of the pseudocritical point to be 

^pscudo ~ L , (39) 

where L = \/S is the characteristic size of the system. Thus from Eqs. (l3Tj) and ( 139]) we find 
that the shift exponent is A = 1 for odd TV. In Fig. 2(a) and 2(b) we plot the t dependence 
of the specific heat for 16 x 17 lattice. We can see from Fig. 2(b) that the position of the 
specific heat peak is shifted from zero. 

For even TV (TV = 2N), one can see from Eqs.®, (1TD]) . (ITT]) . (j25]| and flSO]) that the 
partition function Z M ^{t) and the specific heat C(t, Ai,Af) are an even function with 
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respect to its argument t 

C(t,M,Af) = C(0,M,M) + -C^(0,M,Af) + -C^(0,M,Af) + 0(f). (40) 

The partition function given by Eq. (QQ) is a polynomial function of its argument t. In 
the case of even M. and M the partition function is an even function with respect of its 
argument t and hence the number of diagonal bonds n t also should be even. There is 
also a simple geometrical explanation why, in the case of even-even lattices, the number 
of diagonal dimers should be even. Let us consider an even-even (2Ai x 2J\f) lattice with 
periodic boundary conditions in the vertical and horizontal directions. Such a lattice can be 
divided into two sublattice A and B, as shown in Fig. 3. Each sublattice consists of 2M.M 
sites in a such way that every horizontal or vertical edge connects a site in sublattice A to 
one in sublattice B, while diagonal edges connect two sites within sublattice A or two sites 
within sublattice B. Note that such a division into A and B sublattices is impossible for 
even-odd (2.M x 2M + 1) lattices with periodic boundary conditions, since in that case, one 
can always find a horizontal or vertical edge which connects two sites in the same sublattice 
A or B. Thus, in the even-even case, each horizontal or vertical dimer occupies one site 
from sublattice A and another site from sublattice B, while diagonal dimer occupies two 
sites from sublattice A or B. If one diagonal dimer occupies two sites from sublattice A one 
should have another diagonal dimer which occupies two sites of sublattice B in order to 
insure that remaining sites can be occupied by horizontal and vertical dimers. Thus, in the 
case of even M and N, only even number of diagonal bonds are allowed. Similar geometrical 
considerations for the dimer model on ( M. x M) lattice with free boundary conditions lead to 
the conclusion that for both even-even and even-odd lattices the number of diagonal bonds 
should be even. 

Thus the first derivative of C(t, A4,Af) vanishes exactly at 

^pseudo 0. (41) 

In Fig. 4(a) and 4(b) we plot the t dependence of the specific heat for a 16 x 16 lattice. We 
can see from Fig. 3(b) that the position of the specific heat peak t pseu d is equal exactly to 
zero. Therefore the maximum of the specific heat (the pseudocritical point t psC udo) always 
occurs at vanishing reduced temperature for any finite M x 2N lattice and coincides with 
the critical point t c at the thermodynamic limit. From Eqs. (}3~Tj) and (141!) we find that the 
shift exponent is A = oo for even M . 

12 
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Figure 2: (a) The behavior of the specific heat on M x N lattices with even Ai = 16 and odd 
N = 17. (b) The behavior of the specific heat on the same lattice for small t. 




Figure 3: Division of lattice on two sublattice A (black circles) and B (white circles). 

Thus we have found that the shift exponent A for the specific heat depend in a crucial 
way on the parity of the number of lattice sites along the lattice axis M . For odd M we 
obtain for the shift exponent A = 1, while for even M we have found that the shift exponent 
A is infinity ( A = oo). 



V. CONCLUSION 



We analyze the partition function of the dimer model on A4 xj\f triangular lattice wrapped 
on torus obtained by Fendley, Moessner and Sondhi 22j. From a finite-size analysis we have 
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Figure 4: (a) The behavior of the specific heat onMxJV lattices with even M = 16 and N = 16. 
(b) The behavior of the specific heat on the same lattice for small t. 

found that the dimer model on the triangular lattice can be described by conformal field 
theory having central charge c = —2. Thus we have shown that the dimer model on the 
triangular lattice belongs to the same universality class as the dimer model on the square 
lattice, while the dimer model on the honeycomb lattice belongs to another c = 1 universality 
class. In addition, we have found that the shift exponent A depends in a crucial way on the 
parity of the number of lattice sites along the lattice axis M\ for odd M we obtain A = 1, 
while for even M we have found that A = oo. In the former case, therefore, the finite- 
size specific- heat pseudocritical point is size dependent, while in the latter case it coincides 
with the critical point of the thermodynamic limit. This adds to the catalog of anomalous 
circumstances where the shift exponent is not coincident with the correlation-length critical 
exponent. The present circumstance manifests the additional feature that the shift-exponent 
is boundary-condition dependent. 
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